Abstract 



We have performed a search for the leptonic B decay f^^i^fi with data 

cohected at the T(4S') resonance by the BABAR experiment at the PEP- 
II storage ring. With an integrated luminosity of approximately 81.4 fb^^ 
(88.4 million BB pairs), we find no convincing evidence for a signal and set 
a preliminary upper limit on the branching fraction of B{B~^ — > ^'^v^) < 
6.6 X 10~^ at the 90% confidence level. 
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1 INTRODUCTION 



The purely leptonic decays — > £~^ui { £ = e, fi, or r) proceed through the annihilation (Fig. 
of the quark-antiquark pair in the meson to form a virtual W boson. In the Standard Model(SM), 
the branching fraction can be calculated as ^, 

BiB^ ^ l^u,) = (l - (1) 

where Gi? is the Fermi coupling constant, me and tjib are the lepton and B meson masses, and 
tb is the -B"*" lifetime. The decay rate is sensitive to both the Cabibbo-Kobayashi-Maskawa matrix 
element Vub and the B decay constant /b which describes the overlap of the quark wave functions 
within the meson. Once \Vub\ has been more precisely measured in B semi-leptonic decays, Jb could 
be extracted from a measurement of the B~^ — > i'^i'e branching fraction. Currently, the uncertainty 
on Jb is one of the main factors limiting the determination of \Vtd\ from precision Bm^ mixmg 
measurements. 

The Standard Model estimate of the branching fraction for B~^ — > t^Vt is about 9 x 10^^ 
assuming tb = 1.674 ps, \Vub\ = 0.0036 and Jb = 198 MeV jT]. Due to helicity suppression, the 
expected branching fraction for B~^ — > n'^v^ is reduced to roughly 4 x 10"''. However, physics 
beyond the Standard Model could substantially increase these predictions. Charged Higgs boson 
effects may greatly enhance the branching fraction in certain two Higgs doublet models [21 . Sim- 
ilarly, this decay may be enhanced through mediation by leptoquarks in the Pati-Salam model of 
quark-lepton unification |^. 




Figure 1: The SM B'^ f-'^i^t annihilation diagram. 

Purely leptonic B decays have not yet been observed experimentally. The CLEO [l] and Belle 
Collaborations have set 90% confidence level upper limits on the B~^ — > ^'^v^ branching fraction. 
The most stringent limit is currently from Belle, 

^ li-^v^) < 6.8 X 10"*^ 

usmg 60 fb"^ of data collected at the r{4S) resonance. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The sample consists of an integrated luminosity of 81.4 fb~^ accumulated at the T{4S) 

^Charge-conjugation is implied throughout this paper. 
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resonance and 9.6 fb~ accumulated at a center-of-mass (CM) energy about 40 MeV below the 
T{4:S) resonance. The on-resonance sample corresponds to about 88.4 million BB pairs. The 
collider is operated with asymmetric beam energies, producing a boost of the T(4S) along the 
collision axis of Pj = 0.55 in the laboratory frame. 

The BABAR detector is optimized for the asymmetric-energy beams at PEP-II and is described 
in detail elsewhere Charged particle momentum and direction are measured with a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) which are contained 
in the 1.5 T magnetic field of a superconducting solenoid. Located just outside the DCH, a detector 
of internally reflected Cherenkov radiation (DIRC) provides separation of and vr"'". The energies 
of neutral particles are measured by the electromagnetic calorimeter (EMC) which is constructed of 
6580 CsI(Tl) crystals. The flux return of the solenoid is instrumented with resistive plate chambers 
(IFR) for the identification of muons and K^. Averaged over the momentum and polar angle 
distributions of muons from B~^ fJ-^i^^, the muon identification efficiency is about 61% with a 
pion misidentification probability of about 2%. 

A Geant4 based Monte Carlo (MC) simulation was used to optimize the signal selection criteria. 
A sample of 52,000 simulated B~^B~ events where B^ l^-^^n ^^nd the B~ decayed generically has 
been studied. Background sources considered include e^e~ BB, e^e~ qq {q = u, d, s, and ), 
and e+e" — > r+r" in quantities comparable to the data luminosity. 



3 ANALYSIS METHOD 

B^ — > fi^i^fj, is a two-body decay so the muon must be mono-energetic in the B rest frame. The 
momentum of the muon is given by 

^ ml -ml ^ TUB 
^ 2mB ^ 2 ■ ^ ' 

At BABAR, the CM frame is a good approximation to the B rest frame so we initially select 
well-identified muon candidates with momentum pcM between 2.25 and 2.95 GeV/c in the CM 
frame. Since the neutrino produced in the signal decay is not detected, any other charged tracks 
or neutral energy in a signal event must have been produced by the decay of the companion B. 
Therefore, the companion B can be reconstructed from the remaining visible energy in the event. 
Signal decays can then be selected using the kinematic variables IS.E and energy-substituted mass, 
rriES) defined by 

A-E = Eb — -Ebeam, (3) 

and 

"iES = V^beam " ^ (4) 

where pB and Eb are the momentum and energy of the reconstructed companion B candidate 
in the CM frame and -Ebeam is the beam energy in the CM frame. The reconstruction of the 
companion B includes all charged tracks whose distance of closest approach to the beam-spot is 
less than 1.5 cm in the xy plane and within 10 cm along the beam axis. We also include all 
neutral calorimeter clusters with cluster energy greater than 30 MeV. Particle identification is 
applied to the charged tracks to identify electrons, muons, kaons and protons in order to apply 
the appropriate mass hypothesis to each track and thus improve the A£' and rriES resolution. 
In addition, events with additional identified leptons are discarded to discriminate against events 
containing additional neutrinos. The MC indicates that this requirement removes about 24% of the 
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Figure 2: The distributions of AE and ttt-es for on-peak data and MC after muon candidate 
selection. The signal distributions are overlaid (dashed histograms) with an arbitrary normalization. 



signal decays. Figure |^ shows the distributions of AE and ttt-es for the on-peak data, background 
MC and signal MC after muon candidate selection. For a signal decay, we expect the energy of 
the companion B to be consistent with the beam energy in the CM frame so that AE peaks near 
0. Due to energy losses from detector acceptance, neutral hadrons and additional neutrinos, this 
distribution is shifted toward negative AE. We expect that tti-es should peak near the B mass for 
signal decays. 

Once the companion B is reconstructed, we make a refined estimate of the muon momentum 
in the B "rest" frame. We use the momentum direction of the companion B and assume a total 
momentum of 320 MeV/c in the CM frame (from the decay of the T{4S) BB) to boost the muon 
candidate into the reconstructed B rest frame. Figure El shows the muon candidate momentum 
distribution in the B rest frame, p* , for all muon candidates in the signal MC. The dashed curve 
is the momentum distribution of the same events in the CM frame. 

Backgrounds may arise from any process producing charged tracks in the momentum range of 
the signal, particularly if the charged tracks are true muons. The two most significant backgrounds 
are the B semi-leptonic decays involving b — > ufiv transitions where the endpoint of the muon 
spectrum approaches that of the signal, and non-resonant qq (continuum) events where a charged 
pion is mistakenly identified as a muon. In the continuum events, there must also be significant 
missing energy due to detector acceptance, neutral hadrons, or additional neutrinos that mimic 
the signature of the expected neutrino. We are able to reduce these backgrounds by tightening the 
selection on the muon candidate momentum to 2.58 < p* < 2.78 GeV/c. The cut is asymmetric 
about the signal peak due to the decreasing momentum distribution of the backgrounds. 

Continuum backgrounds are further suppressed using event shape variables. The light-quark 
events tend to produce a jet-like event topology as opposed to BB events which tend to be more 
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Figure 3: The muon candidate momentum distribution in the reconstructed B rest frame for ah 
muon candidates in the signal MC. The dashed curve is the momentum distribution of the same 
events in the CM frame. The arrows indicate the selected signal region. 



spherical. We define a variable, 6 , which is the angle between the muon candidate momentum and 
the thrust axis of the rest of the event in the CM frame. For continuum background, | cos 9'^\ peaks 
sharply near one while the distribution is nearly flat for signal decays. By requiring | cos < 0.55, 
we are able to remove about 98% of the continuum background while retaining about 53% of the 
signal decays. We also use the direction of the missing momentum in the laboratory frame to 
discriminate against continuum backgrounds. The missing 4- momentum is calculated as 

^ pr(4S) _pB _pf,^ 

where p'^(^^) is determined from the beam energies, and and P'^ represent the reconstructed 
companion B and signal muon candidate respectively. In continuum decays, the missing momentum 
is often due to undetected particles that were outside the detector acceptance so that the polar 
angle of the missing momentum, 9'^, lies near the beam axis. Therefore, we require | cos 0*^1 < 
0.88 so that the missing momentum is directed into the detector's fiducial volume. In addition, 
we require that the event contains at least four charged tracks and that the normalized second 
Fox- Wolfram moment W is less than 0.98. Figure |3] shows the on-peak data and MC distributions 
of I cos9^\ and | cos^'^l. The events in these plots have passed the requirement 2.58 < p* < 2.78 
GeV/c. For comparison, the signal MC is overlaid with an arbitrary normalization. The remaining 
BB background is not visible due to the earlier requirement on p* . 

We select B~^ — > /u^t'^ signal events with simultaneous requirements on AE and mES, thus 
forming a "signal box." The dimensions of the signal box, as well as the above requirements on 
p*, I cos 0-^1 and | cos9'^\, were determined using an optimization procedure that searches the cut- 
parameter space to find the combination of cuts that maximizes the quantity S/\/ S + B where S 
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Figure 4: The distributions of | cos^'^l and | cos 0*^1 for on-peak data and MC. The events in these 
plots have passed the requirement 2.58 < < 2.78 GeV/c. The signal distributions are overlaid 
(dashed histograms) with an arbitrary normalization. 



and B are the expected signal and background yields in the MC respectively. The signal branching 
fraction was set to 3x10"'' during the optimization procedure. The resulting box is defined by 
-0.75 < /^E < 0.5 GeV and rriES > 5.27 GeV/c^. The MC suggests that about 25% of signal 
decays passing all previous cuts fall within the signal box. After applying all selection criteria, the 
H^i^fj, efficiency is determined from the MC to be 2.24 it 0.07 % where the uncertainty is 
due to MC statistics. Figure El^a) shows the distribution of AE vs tties for the signal MC. The 
signal box is represented by the solid lines. 

In addition to the signal box, we have defined a slightly larger blinding box. The data within 
the blinding box was kept hidden until the analysis was completed in order to remove any possible 
experimenter's bias. Finally, we define several sideband regions that will be denoted the fit sideband 
and the AE sidebands. The boundaries of these regions in the {AE, uies) plane are listed in tabled 

Table 1: The boundaries of the signal box and various sidebands defined for this analysis. 



region 


range ( GeV) 


rriES range (GeV/c^) 


signal box(SB) 


[ -0.75, 0.50 ] 


> 5.27 


blinding box 


[ -1.30, 0.70 ] 


> 5.24 


fit sideband 


[ -0.75, 0.50 ] 


[ 5.10, 5.24 ] 


AE sideband (bottom) 


[ -3.00, -1.30 ] 


> 5.10 


AE sideband (top) 


[ 0.70, 1.50 ] 


> 5.10 
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We estimate the background in the signal box assuming that the rriEs distribution is described 
by the Argus function JHI . This assumption is consistent with the observed distributions in the MC 
and data AE sidebands. Due to the correlation between AE and tit-es introduced by the inclusive 
reconstruction of the companion B, it is not possible to determine the Argus shape from the AE 
sidebands. Therefore, the Argus parameter is determined from an unbinned maximum likelihood 
fit using the data in the region defined by -0.75 < AE < 0.5 GeV and 5.10 < uies < 5.24 GeV/c^. 
The Argus shape (A) is extrapolated through the signal box and constrained to be at the endpoint 
which is fixed at 5.29 GeV/c^. Figure El shows the results of the fit. The expected background is 
given by, 

^Vbkg — ^Vfit X g 24 . / n = ^^fit ^ ^Argus (.Dj 

J5.10 A{mEs)amES 

where Nf^t is the number of events contributing to the fit. The ratio i?Argus is the area of the 
background probability density function in the signal box divided by the area within the fit region. 
The result is A^bkg = 5.OII4 events. The uncertainty is determined by varying the Argus parameter 
by the itlcr uncertainty from the fit. This technique may underestimate the background since any 
peaking component within the blinding box would not be accounted for. Therefore, any upper 
limit obtained on the branching fraction should be conservative. When this procedure is applied 
to the MC, the resulting background estimate is 5.2 ± 0.5 events, in agreement with the true value 
of 5.7 zb 0.5 events. The MC indicates the background in the signal box after all cuts are applied 
is composed of about 58% qq (where q = u, d, or s), 22% cc and 20% BB events. 



4 SYSTEMATIC STUDIES 

To set an upper limit on the B^ /^^i^^t branching fraction we must evaluate systematic uncer- 
tainties in the luminosity (number of B^ in the sample), the background estimate, and the signal 
efficiency. The number of B^ mesons in the on-peak data sample is estimated to be 88.4 million 
with an uncertainty of 1.1%. The uncertainty in the background estimate is determined by varying 
the Argus shape within the itlo" uncertainty from the fit. 

The uncertainty in the signal efficiency includes the muon candidate selection (particle identi- 
fication and tracking efficiency) as well as the reconstruction efficiency of the companion B. The 
muon identification efficiency has been studied using muon control samples taken from e~^e~ 
e~^e~ events in the data. For muons in the momentum and polar angle region of our signal 

candidates, we estimate the uncertainty in the identification efficiency from the control sample to 
be 4.2%. The tracking efficiency for the muon candidate was evaluated from the fraction of tracks 
reconstructed in the SVT that are also found in the DCH. We find that the tracking efficiency is 
overestimated in the MC by 0.8%. Therefore, we reduce the MC signal efficiency by this amount 
and assume a systematic error of 2%. 

The companion B reconstruction has been studied using a control sample of B~^ — > D^*')^7r~^ 
events. This is also a two-body decay so it is topologically very similar to our signal. Once 
reconstructed, the pion can be treated as if it were the signal muon and the D^*^'^ decay products 
can be ignored to simulate the neutrino. Then the companion B is reconstructed in the control 
sample as it would be for signal. We then compare the efficiencies for each of our companion B 
selection cuts in the B~^ — > D^*^^7r~^ data and MC to quantify any data/MC discrepancies that 
may affect the signal efficiency. We find that the efficiency after all selection cuts is lower in the 
data by a factor of 0.94zb0.04 where the uncertainty is due to the statistics of the data and MC 
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Figure 5: Results of the Argus fit for the background estimate. The fit is performed only on the 
region 5.10 < uies < 5.24 GeV/c^ (diagonally shaded area) and extrapolated into the signal region 
represented by the crosshatched area. 

control samples. The signal efficiency obtained from the MC is therefore corrected by this factor. A 
summary of the systematic uncertainties in the signal efficiency is given in table IS After applying 
the necessary corrections to the MC, we estimate the signal efficiency to be 2.09 it 0.06 (stat) it 
0.13 (syst) %. 

Table 2: Contributions to the systematic uncertainty on the signal efficiency. 



source 


correction 


uncertainty 


tracking efficiency 


0.992 


2.0% 


muon identification 




4.2% 


companion B reconstruction 


0.94 


4.3% 


total 


0.932 


6.3% 



5 PHYSICS RESULTS 

In the on-resonance data we find 11 events in the signal box where 5.0l}'4 background events are 
expected. The probability of a background fluctuation yielding the observed number of events or 
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Figure 6: The distributions of AE vs m-ES in the — > ^^t'^ signal MC (a) and the on-resonance 
data (b). The signal box is represented by the sohd hne while the grand sideband is represented 
by the dashed line. 



more is about 4%. The distribution of the data in the {AE, niEs) plane is shown in figure l^lb). 
Figure 13 shows the p* distribution of the on-resonance data after all other selection cuts have been 
applied. The region between the two arrows (2.58 < p* < 2.78 GeV/c) indicates the selected signal 
candidates. We set an upper limit on the — > /i^t'^ branching fraction using B{B~^ A^^^^^i) < 
uul/S where nuL is the 90% CL upper limit on the number of signal events observed and S is 
the sensitivity of the experiment which is the product of the signal efficiency and the number of 
charged B mesons in the sample. To determine the number of charged B mesons we assume equal 
production of B^ and B~^ in T(4S') decays. Systematic uncertainties are included in the limit 
following the prescription given in reference IH]. The preliminary result is 

B{B+ ii+v^) < 6.6 X 10"^ 

at the 90% confidence level. 
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